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Abstract 

We report preliminary measurements of absolute hadronic branching fractions of D s mesons 
determined using a double tag technique. These measurements are from 195 pb _1 of e + e~ collisions 
recorded at center of mass energies near 4.17 GeV with the CLEO-c detector at CESR. We obtain 
absolute branching fractions for Df decays to KgK + , K~ K + tt + , K~ K + ir + ir° , 7r + 7r + 7r _ , ir + g, 
and Tr + n'. We discuss the problems inherent in measuring accurately the branching fraction for 
Df — > (f>ir + , which is often used as a reference mode for measurement of other Df branching 
fractions, and provide a measurement of a branching fraction that may be useful for this purpose. 



* Submitted to the 33 rd International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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I. INTRODUCTION 



Absolute measurements of hadronic charm meson branching fractions play a central role 
in the study of the weak interaction because they serve to normalize many important B 
and D meson branching fractions. The importance of many of these branching fractions is 
enhanced by their roles in determining elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [1]. 

Absolute D s l branching fractions have not been measured as accurately as those for D + 
and D° [2], due to difficulties in determining the total number of D s mesons produced in 
an experiment. The current best measurements of B(Df — > 07r + ), which is usually used to 
normalize measurements of other D s decay modes, are CLEO and BaBar results obtained 
using B° -> D*-D* + events and a missing mass technique to extract the number of Df 
[3,4]. 

In this contribution we describe measurements of absolute D s branching fractions using 
a double tag technique. This technique was pioneered by the MARK-III Collaboration [5, 6] 
and more recently exploited in our precision measurements of absolute D + and D° branching 
fractions [7]. This technique relies on the fact that charm or strange quarks are always 
produced with the corresponding antiquark (cc or ss) in e + e~ annihilation. Hence, below 
the thresholds for DfD~K° and DfD°K~ near E cm = 4.33 GeV, production of a Df is 
always associated with production of a Dj , so observation of a Df can serve as a tag for 
production of a Df. In this experiment, data were accumulated near E cm = 4.17 GeV (below 
the threshold for Df^Df^ production) where the production of DfD~ pairs is dominated 
by e + e~ — > Df^Df [8] with the subsequent decays D* ± — > Df^j or Dfn°. 

Following the MARK III Collaboration, events in which either a Df or Df is recon- 
structed without reference to the other particle are called single tag (ST) events, while 
events in which both the and Df are reconstructed are called double tag (DT) events. 
Absolute branching fractions can then be determined from the ratios of DT events to ST 
events. If CP violation is negligible, then the branching fractions Bj and Bj for Df — > j and 
Df — > j are equal. However, the efficiencies q and ej for detection of these modes may be 
somewhat different since the cross sections for scattering of particles on nuclei may depend 
on the signs of the particles. With the assumption that Bj = Bj, the numbers JVj and Nj of 
Df — > % and Df — > j ST events will be, 

Ni = N D * Ds Bi€i and Nj = N D * Ds Bj€j, 

where Nd*d s is the number of Df^Df events produced in the experiment. Then, the number 
of DT events with Df — > i and D~ — > j will be, 

N ij = -V/r • 

where is the efficiency for detecting the DT. Hence, the ratio of DT events (N^) to ST 
events (Nj) provides an absolute measurement of the branching fraction £>j, 

N-- e-- 



Generally D s will refer to either Df or D s , and specification of an explicit Df state and its decay 
daughters will imply a corresponding relationship for the Df and its daughters. 
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Note that eq pa e^ej, so many systematic uncertainties nearly cancel in the ratio ej/e^j, and a 
branching fraction Bi obtained using this procedure is nearly independent of the efficiency of 
the tagging mode. Of course, Bi is sensitive to e« and its uncertainties. Finally, the number 
of D+D~ pairs that were produced is given by, 

N-N- e - 

"did. - -j^r — • 

We studied six decay modes, D+ -> K° S K + , D+ -> K-R + tt + , D+ -> fT-fT+Tr+Tr , 
/}+ — > 7r + 7r + 7r _ , — -> 7i + i], and -D+ — ► and their charge conjugates in this analysis. 

For each mode we report the average of the branching fractions for and D~ . We detected 
neutral kaon candidates in the K s — > 7r + 7r~ mode, i] candidates in the r\ — > 77 mode, and 
7/ candidates in the 7/ — > 7r + 7r~77 mode. 



II. THE CLEO DETECTOR AND PARTICLE SELECTION 

The CLEO-c detector is a modification of the CLEO III detector [9-11], in which the 
silicon-strip vertex detector was replaced with a six-layer vertex drift chamber, whose wires 
are all at small stereo angles to the chamber axis [12]. The charged particle tracking system, 
consisting of this vertex drift chamber and a 47-layer central drift chamber [10], operates 
in a 1.0 T magnetic field, oriented along the beam axis. The momentum resolution of 
the tracking system is approximately 0.6% at p — 1 GeV/c. Photons are detected in an 
electromagnetic calorimeter, composed of 7800 CsI(Tl) crystals [9], which attains a photon 
energy resolution of 2.2% at E 7 — 1 GeV and 5% at 100 MeV. The solid angle coverage 
for charged and neutral particles of the CLEO-c detector is 93% of 4tt. We utilize two 
particle identification (PID) devices to separate from n ± : the central drift chamber, 
which provides measurements of ionization energy loss (dE/dx), and, surrounding this drift 
chamber, a cylindrical ring-imaging Cherenkov (RICH) detector [11], whose active solid 
angle is 80% of 4tt. The combined dE/dx-KlCH PID system has a pion or kaon efficiency 
> 90% and a probability of pions faking kaons (or vice versa) < 5%. 

The response of the detector is studied with a detailed GEANT-based [13] Monte Carlo 
(MC) simulation of the CLEO detector for particle trajectories generated by EvtGen [14] 
and final state radiation (FSR) predicted by PHOTOS [15]. Simulated events are processed 
in a fashion similar to data. The integrated luminosity (£) of the data sample is measured 
using e + e~ Bhabha events in the calorimeter [16], where the event count normalization is 
provided by the detector simulation. 

Charged tracks are required to be well-measured and to satisfy criteria based on the track 
fit quality. They must also be consistent with coming from the interaction point in three 
dimensions. Pions and kaons are identified by consistency with the expected dE/dx and 
RICH information, when available. We select K s candidates from pairs of oppositely-charged 
tracks that are constrained to a vertex, and that have \M(ti + it~) — M K o\ < 6.3 MeV/c 2 , 
where M K o [2] is the known K s mass. The 6.3 MeV/c 2 limit is approximately 3a. 

We form 7r° (rj) candidates from photon pairs with invariant mass within 3 standard 
deviations (a) of the known 7r° (77) mass, M n o (M„) [2], where a m 5-7 MeV/c 2 depending 
on photon energy and location. These candidates are then fit kinematically with their masses 
constrained to M w o (M v ). We reconstruct rj' candidates in the decay mode 7/ — > 7t + 7t~t7 and 
require 0.9478 < M(rf) < 0.9678 GeV/c 2 . 
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III. DATA SAMPLE AND CANDIDATE SELECTION 



The data sample for this measurement consists of 179 pb 1 collected at E cm = 4.17 GeV, 
10 pb" 1 collected at E cm = 4.16 GeV, and 6 pb" 1 collected at E cm = 4.18 GeV. We 
collected the 4.16 GeV data and 4.18 GeV data in a scan of center-of-mass energies designed 
to determine the optimal energy for studying D s decays [8]. 

Most of the D s mesons at E cm near 4.17 GeV come from D^DJ events, and we impose 
kinematic cuts on the D s candidates to require this topology. In such events, D*^ — > Dfj 
or Dfn , so one of the D s is produced directly and the other is a secondary from D* decay. 
We use two variables to select D s decay candidates and isolate signals from backgrounds. 
One variable is the candidate invariant mass, 



where p(D s ) and E{D S ) are the total momentum and energy, respectively, of the daughters 
of the candidate D s . The particle energies are computed from particle momenta using the 
masses determined by PID. The other variable is the beam constrained mass, 



obtained by substituting the beam energy E for the energy E(D S ) of the candidate. We 
make the following observations and choices based on these observations for this analysis: 

• We do not require full reconstruction of the event; i.e., we make no attempt to find 
the 7 or ir° resulting from the D* decay. This avoids a large efficiency loss and avoids 
substantial contributions to systematic errors that would arise from soft photon detection 
and which would affect ST and DT yields differently. 

• The value of E cm determines the momenta of the D s and D*, so the momentum of the 
direct D s is tightly constrained, and the momentum of the D s from D* decay has a 
somewhat wider distribution about the same value. 

• The M bc distribution of direct D s mesons has a peak about 0.010 GeV/c 2 wide centered 
at 2.04 GeV/c 2 - somewhat above M Ds = 1.9683 GeV/c 2 [2]. The M bc distribution of 
the D s daughters from D* decay is broader and relatively flat from about 2.015 GeV/c 2 
to about 2.065 GeV/c 2 . These two components of the M bc mass distribution are clearly 
evident in Fig. 1. The upper limit of the M bc distribution is just the beam energy E , 
while the lower limit is determined by the maximum momentum of the secondary D s 
candidates. 

• In our selection of D s candidates, we use M bc as a proxy for this momentum, since it 
is insensitive to variations in E cm and allows easy integration of data taken at different 
energies. 

• We use the invariant mass M(D S ) of the D s candidate as our primary analysis variable. 
It has very little correlation with M bc in the region of interest. 

We require \M(D S ) - M Ds \ < 0.085 GeV/c 2 and a "loose" cut M bc > 2.01 GeV/c 2 , 
for D s candidates in all decay modes. For the D+ — > K~ K + -K + ir° and Df — > it + it + it~ 
modes, we require additional "tight" cuts of M bc > 2.034 GeV/c 2 and M bc > 2.035 GeV/c 2 , 
respectively. The loose cut includes the entire range of M bc resulting in maximum efficiency. 
The tight cut selects all of the directly-produced D s and over half of the indirect D s as well. 
This reduces the background while reducing the impact of initial-state radiation smearing 
of efficiency. 




(1) 




(2) 
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Furthermore, for ST candidates we require: 

• Df — > K S K + : We reject candidates with M(K + n~) > 1.83 GeV/c 2 to eliminate a small 
background from D° — > i^ + 7r~ decays. 

• L>+ -> K~K + n + : We reject candidates with 1.845 < M(K~K + ) < 1.88 GeV/c 2 to 
eliminate a small background from D° — > K~K + decays. 

• L>+ — > K~K + it + it : We require p^o > 0.1 GeV/c to reduce a background from D*^ 
and £>*° decays, and we reject candidates with 1.86 < M(K~ K + n + ) < 1.88 GeV/c 2 to 
eliminate a small background from D + — > i^~fT + 7r + decays. 

• — > 7T + 7r + 7r _ : We reject candidates if either 7r + 7r~ candidate has 475 < M(tt + tt~) < 
520 MeV/c 2 (to eliminate K s — > 7r + 7r~ decays), or either candidate has 1.840 < 
M(ir + ir~) < 1.885 GeV/c 2 (to eliminate D° — > 7r + 7r~ decays). Also, when the 7r~ is 
assumed to be a K~ and is assigned the K~ mass, we reject the candidate if either 
K~n + combination had 1.845 < M(K~tt + ) < 1.880 GeV/c 2 to eliminate background 
from D° -> K^7T+. 

• In every event, for each mode and each charge, we choose the best ST candidate. We 
accept the candidate with the smallest value of \M bc — M|, where M = 2.040553 GeV/c 2 
is the value at which M bc peaks for daughters of direct D s decays. An event can have 
candidates in more than one mode for a given D s charge. 

Furthermore, for DT candidates we require: 

• In a DT event, one of the D s will be direct and the other will be the daughter of a D* s 
decay. We impose the loose M bc requirement on both of the D s candidates. 

• Each of the D s candidates in a DT event must pass the ST requirements described above. 

• To resolve multiple candidates, we select the candidate with total invariant mass closest 
to 2My) s . This variable is orthogonal to the difference between the masses of the two 
candidates, which is used in selecting signal and background sideband events as described 
below. 

IV. EFFICIENCIES AND YIELDS 

To develop our analysis procedures and estimate efficiencies, we simulated two types of 
D s meson decays: 

• "signal" decays in which a or D~ decays in one of the 6 modes measured in this 
analysis, and 

• "generic" decays in which a Df or D~ decays in a manner consistent with the decay 
modes and branching fractions in the PDG compilation [2], extended with additional 
decays that are expected to be present. 

Using these two types of simulated decays, we generated three types of Monte Carlo (MC) 
events: 

• generic Monte Carlo events, in which both the Df and the D~ decay generically, 

• single tag signal Monte Carlo events, in which either the or the D~ always decays in 
one of the 6 modes measured in this analysis while the D~ or Df, respectively, decays 
generically, and 

• double tag signal Monte Carlo events, in which both the D+ and the D~ decay in one of 
the 6 modes in this analysis. 
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In these events, both the Df and D~ have a 50% probability of being produced from a Df*. 

We also generated MC samples of continuum u, d, s quark production and production of 
non-strange charmed meson states, which form the major background to this analysis. The 
relative production rates are set by information obtained from the scan of this energy region 
[8]- 

We obtained the ST efficiencies and yields by fitting M(D S ) distribution for each mode 
in single tag MC events. The fit functions and their parameters, as well as efficiencies were 
determined from M(D S ) distributions from signal Monte Carlo simulations. 

The ST distributions are fit with a linear backgrounds and a signal lineshapes. Two 
Gaussian distributions were chosen for the lineshapes for the K$K+, K~K+tt+, 7r + 7r + 7r~, 
and tt+i]' modes. The lineshapes chosen for K~ K + 7t + tt and n+i] decays was the sum of 
a narrow Gaussian and a wider Crystal Ball function [17], whose low mass tail simulates 
energy leakage in the Csl calorimeter. The lineshape parameters are fixed separately for 
each charge from candidates that are matched to the generator- level D s . The ST efficiencies 
obtained from signal MC events are given in Table I. The fits to M(D S ) distributions in data 
are illustrated in Figs. 2-7, and the corresponding yields from data are given in Table II. 



TABLE I: Single tag efficiencies (in %) obtained from signal MC events. Sec. VI describes small 
corrections that must be applied to these efficiencies. 







Df Mode 


Df 


Ds 






K S K+ 


37.58 ± 0.24 


37.43 ± 0.24 






K~K+tt+ 


44.79 ± 0.25 


43.90 ± 0.25 






K-K+TT+TT 


12.32 ± 0.19 


12.69 ± 0.19 






tt+tt+tt^ 


49.52 ± 0.26 


49.36 ± 0.26 






tt+t] 


19.99 ± 0.21 


19.09 ± 0.21 






tt+t]' 


5.29 ± 0.12 


5.53 ± 0.12 




TABLE II: S 


ingle tag yields obtained from fits 


to the M(D S ) distributions of data. 




K S K+ 


K-K+-K+ 




TT+TT+TT^ TT+T] TT+7]' 


Df 


1055.0 ± 39.4 


4319.1 ± 88.8 


1342.3 ± 101.9 


969.5 ± 79.3 547.0 ± 49.6 361.8 ± 23.4 


Ds 


927.6 ± 37.3 


4352.2 ± 88.7 


1371.9 ± 102.9 


946.7 ± 77.6 570.0 ± 49.8 371.5 ± 23.6 



We obtain double tag yields by counting the numbers of events in signal regions in 
M{D~) vs M(Df) planes and subtracting the numbers of events in sideband regions of the 
corresponding planes. We find that combinatoric backgrounds from D+, D°, and continuum 
events tend to be distributed along lines of constant total mass, Mt = M(Df) + M(Dj). 
Both signal and sideband regions use subsets of the My region, 3.915 < Mt < 3.960 GeV/c 2 . 
Within this Mt region, the signal region is \ AM\ < 30 MeV/c 2 , while the sideband regions 
are 50 < \AM\ < 150 MeV/c 2 , where AM = M(Df ) - M(Dj). We see no evidence for 
any background to peak in the signal region, so we are confident that the events in these 
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sideband regions provide reasonable estimates of the backgrounds in the signal regions. DT 
efficiencies obtained from signal MC events are given in Table III and DT yields obtained 
from data are given in Table IV. 



TABLE III: Double tag efficiencies (in %) obtained from signal MC events. Sec. VI describes small 
corrections that must be applied to these efficiencies. 





K S K- 


K+K--K- 


K+K-tt-tt 


TT TT TT + 


7r rj 


7r rj' 


K S K+ 


14.02 ± 0.17 


17.30 ± 0.19 


6.74 ± 0.13 


22.84 ± 0.21 


6.64 ± 0.12 


1.96 ± 0.07 


K-K+K+ 


16.75 ± 0.19 


20.35 ± 0.20 


7.65 ± 0.13 


27.91 ± 0.22 


8.01 ± 0.14 


2.39 ± 0.08 


K-K+TT+TT 


6.74 ± 0.13 


7.88 ± 0.13 


2.40 ± 0.08 


11.09 ± 0.16 


3.17 ± 0.09 


1.00 ± 0.05 


tt + tt + tt~ 


22.48 ± 0.21 


27.69 ± 0.22 


11.22 ± 0.16 


35.94 ± 0.24 


10.91 ± 0.16 


3.17 ± 0.09 


TT + 7] 


6.68 ± 0.12 


8.34 ± 0.14 


3.21 ± 0.09 


10.88 ± 0.16 


3.18 ± 0.09 


0.92 ± 0.05 


TT + 7]' 


1.94 ± 0.07 


2.47 ± 0.08 


0.94 ± 0.05 


3.32 ± 0.09 


0.92 ± 0.05 


0.32 ± 0.03 



V. BRANCHING FRACTIONS 

We are unable to use the standard branching fraction fitter [18] created for our analysis of 
D + and D° branching fractions [7] . This fitter uses x 2 minimization, which is not appropriate 
for the small yields and backgrounds in the DT modes. For this analysis, we have developed 
a binned-likelihood hybrid fitter which utilizes Gaussian statistics for the ST modes and 
Poisson statistics for the DT modes. The likelihood is maximized with the number of D^DJ 
mesons (N D * Da ) and the branching fractions for the individual modes (£>j) averaged over D + 
and Dj as the free parameters. The hybrid fitter treats properly the correlated systematic 
errors among quantities affecting ST efficiencies only. However, unlike the standard fitter, 
this fitter does not handle correlated systematic errors in efficiencies among double tag modes 
or between single tag modes and double tag modes. On the other hand, the approximate 
cancellation of systematic errors between DT and ST efficiencies is still utilized in computing 
overall systematic errors. 

TABLE IV: Double tag yields obtained from data by subtracting backgrounds estimated from 
sidebands. 





K S K- 


K + K~tt- 


K+K-tt-tt 


TT TT 7T + 


TT T] 


tt rj' 


K S K+ 


7.7 


27.0 


18.7 


7.3 


4.0 


5.0 


K-K+U+ 


18.0 


104.7 


43.7 


30.7 


12.0 


8.0 




8.7 


35.7 


14.0 


13.3 


1.0 


5.7 


tt + tt + tt~ 


3.3 


22.7 


16.0 


13.3 


4.7 


4.0 


ir + 7] 


0.0 


10.0 


2.7 


6.0 


1.0 


1.7 


ir + r]' 


3.0 


10.0 


3.0 


3.7 


1.0 


0.0 
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We tested this hybrid fitter extensively in MC simulations and with data from our D + and 
D° analysis. This fitter produces results very close to the standard fitter in that analysis. 
We also tested this hybrid fitter with toy MC samples with ST and DT yields comparable 
to those in this analysis. We find that the hybrid fitter shows less bias and produces pull 
distributions that are more nearly Gaussian and with standard deviations closer to 1. 

The branching fractions and statistical errors resulting from applying the hybrid fitter to 
the ST and DT yields are given in Table V. The table also includes the systematic errors 
described in the next section. 

TABLE V: Branching fractions obtained from the ST and DT yields, using the hybrid fitter. The 
errors are statistical and systematic, respectively. 

Df Mode B (%) 

K S K + 1.50 ±0.09 ±0.05 

K-K+7T+ 5.57 ±0.30 ±0.19 

K-K+ir+ir 5.62 ±0.33 ±0.51 

tt+tt+tt- 1.12 ±0.08 ±0.05 

vr+77 1.47 ±0.12 ±0.14 

vr+r?' 4.02 ±0.27 ±0.30 



VI. SYSTEMATIC ERRORS 
A. Detector Simulation 

Using a missing mass technique, we measure efficiencies for reconstructing tracks, K s de- 
cays, and 7T° decays in both data and simulated events. We fully reconstruct ?/>(3770) — > DD, 
ip(2S) — > J/ipn + 7r~, and tp(2S) — > J/ipir°ir° events, leaving out one particle, for which we 
wish to determine the efficiency. The missing mass of this combination peaks at the mass of 
the omitted particle, whether or not it is detected. Then, the desired efficiency is the frac- 
tion of this peak with this particle correctly reconstructed. We find good agreement between 
efficiencies in data and simulated events. The relative uncertainties in these determinations, 
0.35% per charged particle track, 1.1% per Kg, and 5.0% per n° or rj. Extrapolation from 
the relatively low energies of the pions in ip(2S) — > J/ip^n decay to the higher energies 
encountered in our D s decay modes contributes significantly to the larger uncertainties for 
7T° and rj simulation. 



B. Charged Particle Identification 

We study the simulation of the PID efficiencies for charged pions and kaons using decays 
with unambiguous particle content, such as D° — > K s 7r + n~ and D + — > K^-k + -k + . We find 
small momentum dependent corrections to efficiencies for pions and kaons. The uncertainties 
in these corrections are estimated to be less than 0.2% for pions and 0.3% for kaons. We 
use these corrections to compute efficiency corrections which depend on the decay mode, 
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TABLE VI: Corrections to the efficiencies from signal MC simulations, for particle ID, resonant 
substructure, and their product. The ST efficiencies given in Table I are corrected by the product 
factors in this table, while the DT efficiencies given in Table III are corrected by the product 
factors for the two modes in the DT. The uncertainties are the contributions of these factors to 
the systematic error. The quoted uncertainties for the K~K + ir + and K~ K + tt + 'k® modes cover 
the ranges of corrections for the submodes. 



Mode 


PID 


Resonant 


Product 






Substructure 






0.989 ± 0.003 


1 


0.989 ± 0.003 


K-K+1T+ 


0.971 ± 0.008 


1.000 ±0.015 


0.971 ±0.017 


K-K+1T+1T 


0.961 ±0.014 


1.033 ± 0.060 


0.993 ± 0.062 


7r + 7r + 7r~ 


0.989 ± 0.006 


1.000 ± 0.020 


0.989 ±0.021 


■K + rj 


0.998 ± 0.002 


1 


0.998 ± 0.002 


■K + rj' 


0.968 ± 0.006 


1 


0.968 ± 0.006 



and on submodes for modes with multiple resonant states. The correction factors and 
their uncertainties are listed in Table VI. The quoted uncertainties for the K~K + n + and 
K~K + n + -K° modes cover the range of correction factors estimated for the submodes. 

C. Resonant Substructure 

The resonant substructures of the inclusive multibody modes in this analysis are poorly 
understood. Since final-state momenta and angular distributions are sensitive to the inter- 
mediate resonances in these decays, the efficiencies will depend on the resonant substructure. 
We estimate possible corrections to efficiencies and the associated systematic uncertainties 
by estimating the number of events in our data samples in the different resonant modes. We 
then estimate the range of possible corrections from the corrections for these different reso- 
nant modes. This correction does not apply to the two-body modes K$K + , tt + i], and n + r]'. 
The estimated corrections to the efficiencies and the associated systematic uncertainties for 
the multibody modes are included in Table VI. 

D. Single Tag Fitting Functions 

To study possible differences between the MC lineshape and the lineshape in data, we 
perform two checks: We allow the overall mass and width of the signal lineshape to float 
(keeping the relative normalization and ratio of widths of the two Gaussians fixed), and we 
allow a quadratic term in the background. The maximum excursions of the yields of two 
additional fits from the yields of the normal fits are taken to be the systematic uncertainties. 
These excursions and the systematic uncertainties are included in Table VII. 
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TABLE VII: Systematic uncertainties due to uncertainties in the ST fitting functions. The middle 
two columns are yield excursions (averaged over charge) for using a quadratic background or 
allowing the M{D S ) width and mass to float. The rightmost column is the systematic uncertainty 
assigned to single tag yields in this mode. All numbers are in percent (%). 



Mode 


Second-order 


Float 


Systematic 




Polynomial 


Width & Mass 


Uncertainty 




-0.1 


+0.1 


±0.1 


K-K+TT+ 


-0.5 


+2.0 


±2.0 


K-K+TT+TT 


+8.7 


-10.8 


±10.8 


7r + 7r + 7r~ 


+1.4 


-2.1 


±2.1 


ir + 7] 


-8.8 


+11.1 


±11.1 


ir + 7]' 


-4.8 


+4.5 


±4.8 



E. Initial State Radiation 

Initial state radiation (ISR) can impact the M bc requirement used in selection of decay 
modes by reducing the momenta of the D s candidates, leading to high-side tails on the 
distributions for primary and secondary D s mesons. By definition, M bc is limited by the 
beam energy E 

For the loose Mf, c requirement, we accept all events from the lower kinematic limit to 
the beam energy (i.e. momenta from the upper kinematic limit to zero), so the cut is fully 
efficient regardless of ISR. The tight cuts used for the K + K~it + it and it + it + it~ modes 
eliminate most of the lower half of the Mb c distribution for secondary D s from D* decay. 
Since the effect of ISR on the M^ c distribution is independent of decay mode, we study 
the comparison between MC and data using the high-statistics D+ — > K + K~n + mode. 
First, we compare loose and tight yields in data and MC simulation without ISR. We find 
Roata = Ul/ut = 0.807 ± 0.008, where yi and i/t are the yields in data from the loose 
and tight cuts respectively. The corresponding number from MC simulation without ISR 
is Rjsr = 0.814 ± 0.003, so the two ratios agree within 1%. We then study this ratio in a 
MC simulation with ISR, and find Risk = 0.806 ± 0.005, agreeing with the other two ratios 
within 1%. We do not correct for this effect, but we include a 1% systematic uncertainty in 
the ST efficiencies for the K + K~ir + ir° and it + it + it~ modes where we apply this tight M bc 
cut. 

F. Multiple Candidates 

If there is more than one candidate for a given mode, and we select only one, there is some 
inefficiency associated with making the incorrect choice for this mode. This inefficiency is 
modeled in MC simulations, but if the multiple candidate rates are significantly different in 
data and MC simulations, there is likely to be an efficiency difference. We examine signal and 
sideband multiplicity distributions; the invariant mass regions for these two data samples 
are 1.955 < M(D S ) < 1.980 GeV/c 2 and 1.90 < M(D S ) < 1.93 GeV/c 2 , respectively. We 
compare multiplicity distributions for data and MC simulations after subtracting sideband 
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distributions from signal distributions. There is good agreement between data and MC 
multiplicity distributions, although data multiplicities are generally marginally larger. 

If we assume we can do no better than randomly choose candidates (so, for example, the 
efficiency is down 50% for two candidates, 67% for three), we can limit the effect of this 
difference on the total efficiency. We find that this effect on ST efficiencies is negligible for 
the KgK + , K + K~7i + , and 7i + i]' modes. For the other three modes we assign systematic 
uncertainties to the ST efficiencies of: 3% for K + K~7i + 7r°, 1% for 7r + 7r + 7r~, and 0.5% for 

G. Other Possible Systematic Uncertainties 

Several potential systematic errors do not contribute significantly to the systematic un- 
certainties of these measurements. 

We checked our analysis procedure, the internal consistency of our Monte Carlo simula- 
tions, and the utility of the hybrid fitter by analyzing a generic MC sample corresponding 
to approximately six times the luminosity of the current data sample. This sample included 
appropriately weighted samples of all open charm channels at 4.17 GeV. There is no signifi- 
cant evidence for bias in the analysis or fitting procedure, so we do not include a systematic 
error for our procedures. 

Data were accumulated in three separate data sets; we do not observe significant varia- 
tions of yields per unit luminosity over these data sets. 

The ST yields for and D~ are consistent, so we do not include a contribution to the 
systematic error for the differences in efficiencies or yields for and D~. 

The Df daughters from the isospin violating decay D* + — > D+ir° have a different mo- 
mentum distribution from those from D* + — > -D+7 decay. The yields of the K + K~7t + tt° and 
7r + 7r + 7r~ modes are sensitive to this momentum difference due to the tight Mf, c requirements 
for these modes. However, the effect is small due to the small branching fraction for the 
.D+7r mode. We estimate that the possible effect is of (9(0.14%) and we ignore it. 

Potential differences in M bc resolution in MC simulations and data can affect the efficiency 
for the tight cuts used for the K + K~n + ir° and 7r + 7r + 7r~ modes. The M bc resolution <JM bc 
is dominated by the beam energy spread. If we vary this by 10%, we expect an efficiency 
change of 0(0.2%) or less. 

Because we use invariant mass as our fit variable, the only backgrounds that we expect to 
peak in our signal regions of M(D S ) are those that produce the same final state daughters. 
The Cabibbo suppressed decay Df — > Kg7r + can fake — > 7r + 7r + 7r~ decays. As described 
above, we eliminate this background with cuts on the 7i + tt~ mass for either pair, so we do not 
include a correction or systematic error for this possible background. Via misidentification 
of the K + , the decay Df — > K + it + it~ can fake the KgTc + mode. We search for such a signal 
in data and find none, so we do not include a correction or systematic error for this potential 
background. 

H. Summary of Systematic Errors 

Table VIII summarizes the systematic uncertainties for the inclusive branching fractions 
measured in this analysis. 
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TABLE VIII: Summary of systematic uncertainties in efficiencies for inclusive branching fraction 
measurements. This list excludes systematics from MC statistics and fit function systematics 
obtained in the data fits. Where ranges are given, the values for individual modes can be found 
in Table VI for Particle Identification, Table VII for ST Fit Functions, and Sec. VI F for Multiple 
Candidate Rates. 



Source 


Uncertainty (7oJ 


Affects 


Detector Simulation 


0.35 


tracks 




1.1 


K s 




5.0 


n 




5.0 


V 


Particle Identification 


0.3-1.4 


Per mode 


Resonant Substructure 


1.5 


K+K-TT+ 




6.0 


K-K + 7r+ir 




2.0 


ir + ir + ir~ 


ST Fit Functions 


0.1-11.1 


Per mode 


Initial State Radiation 


1.0 


K-K+IT+TT 




1.0 


7T + 7T + 7r _ 


Multiple Candidate Rates 


0.0-3.0 


Per mode 



The systematic uncertainties discussed in this section were added in quadrature to obtain 
the systematic uncertainties given in Table V for the six inclusive modes included in this 
analysis. 

VII. PARTIAL BRANCHING FRACTIONS FOR Dg K-R+7T+ 

The decay D+ — > 07r + — > K~K + n + , being one of the largest and easiest to reconstruct 
D s decays, has historically been used as a reference mode. Many analyses take advantage of 
the small width and impose cuts on the K~ K + system in a K~ K + n + event to enhance the 
signal significance, generally in the form of an invariant mass cut around the <fi and possibly 
a helicity angle cut (to take advantage of the polarization). Our inclusive branching 
fraction results do not directly connect to analyses that select this submode for background 
suppression. 

We have investigated the K~K + invariant mass spectrum in K~K + tt + decays. The mass 
distribution, obtained by fitting for the D s signal in bins of M(K~ K + ), is shown in Fig. 10. 
Intrinsic detector resolution on M(K~K + ) near the mass is about 1 MeV/c 2 . We observe 
a clear peak arising from the 07r + intermediate state. However there is a broad signal 
from another source, extending from K~K + threshold to masses above M^. This signal 
is found to have kinematic distributions indicative of a scalar. The E687 collaboration 
has published [19] and the FOCUS collaboration has reported [20] Dalitz plot analyses 
of Df — > K~K + n + and find significant scalar contributions from / (980) (or a (980)) in 
this mass region. Because of this extra signal, an exclusive 07r + branching fraction is not 
necessarily useful for normalization purposes, and ignoring the scalar contribution leads to 
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TABLE IX: Efficiency, yields in data, and partial branching fraction results for two choices of 
M(K~K + ) cuts in K~K + ir + events. 



Branching Fraction 


Efficiency (%) 


Data Yield 


Result (%) 


B w 


43.3 ± 1.3 


3062 ± 58 ± 11 


1.98 ± 0.12 ± 0.09 


B 20 


43.4 ± 1.3 


3482 ± 63 ± 13 


2.25 ± 0.13 ± 0.12 



the branching fraction being ill-defined by 0(5%) depending on the specific choice of cuts. 

We report here measurements of the inclusive branching fraction for K~K + tt + decays 
where the K~K + system lies in specific mass ranges. We define 

Bam = B(D+ -> K~K + 7r + ) for \M(K~K + ) - 1.0195 GeV/c 2 | < AM MeV/c 2 

and measure Bw and £>2o- 

To measure these partial branching fractions we perform fits to K~K + n + signal events 
where M(K~ K + ) passes the respective cuts to obtain yields N&m- The efficiencies 
for event reconstruction are obtained from signal MC samples of pure 07r + and /o(980)7r + 
production. The uncertainty arises largely from differences in the reconstruction efficiencies 
for scalar and vector angular distributions. Uncertainties on the fit functions and detector 
resolution smearing of the mass distribution contribute 0.4% to the systematics on the yields 
for both branching fractions. We obtain the partial branching fractions from 

„ N AM 

AM ~ M 

utilizing the value of N D * Da we obtain in the fit to the six inclusive modes. 

The efficiencies, data yields, and partial branching fraction results are shown in Table IX. 

VIII. RESULTS AND CONCLUSIONS 

We have reported preliminary absolute measurements of the branching fractions of D s 
mesons to six final states, with total relative uncertainties ranging from 6.4% to 12.5%. We 
have also obtained measurements of two partial branching fractions, to K~ K + n + substates 
with restricted values of M(K~ K + ). 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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FIG. 1: M hr distribution for D± -> 7T events in data. 
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FIG. 2: Fits to single tag data for Df -> i^g J fsT ± . 
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FIG. 3: Fits to single tag data for Df -» K-R + tt ± . 
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FIG. 4: Fits to single tag data for Df -> if-K+vrV 
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FIG. 5: Fits to single tag data for Df- —>■ ir + ir ir^. 
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FIG. 6: Fits to single tag data for Df — > ir^r]. 
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FIG. 8: Scatter plots of M(D~) vs M(D+) for (left) all six D s decay modes and (right) Df -> 
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FIG. 9: Projections of M(D~) vs M(D+) scatter plots on the AM = M{Df) - M(D~) axis for 
(left) all six D s decay modes and (right) Df — > K~K + tt ± . 
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FIG. 10: M{K + R-) distribution for D+ -> events with M(K + K~) near the mass of 

the 6. 
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